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We report exact calculations of magnetic and superconducting pair-pair correlations for the half- 
filled band Hubbard model on an anisotropic triangular lattice. Our results for the magnetic phases 
are similar to those obtained with other techniques. The superconducting pair-pair correlations at 
distances beyond nearest neighbor decrease monotonically with increasing Hubbard interaction U 
for all anisotropy, indicating the absence of frustration-driven superconductivity within the model. 
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The possibility of superconductivity (SC) in the doped 
square-lattice Hubbard model continues to be discussed 
in the context of cuprate superconductors. Theoretical 
results for the strength of superconducting pair correla- 
tions in this model are controversial, and arguments for 
P, Q and against [1, |3| long range order (LRO) both 
exist in the literature. Yet another class of supercon- 
ductors which have been discussed using the Hubbard 
model are the organic charge-transfer solids (CTS) [H, 
which are known to have strong electron-electron inter- 
actions @. SC in the CTS is reached from proximate 
insulating states not by doping, but by application of 
pressure at constant carrier density. It is therefore natu- 
ral to expect that a simple change of parameters within 
a model appropriate for the insulating states of the CTS 
can give SC. The proponents of spin fluctuation and res- 
onating valence bond models for SC in the CTS sug- 
gest that precisely such an insulator-to-superconductor 
transition can occur within the ^-filled band triangu- 
lar lattice Hubbard model with increasing frustration 
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The frustrated Hubbard model has been applied most 



commonly to the CTS belonging to the the k-(BEDT- 
TTF)2X (k-(ET)2X) family with the highest critical tem- 
perature. In these materials, the ET molecules form 
anisotropic triangular lattices of dimcr pairs 3. While 
K-(ET)2X is formally |-electron filled (|-hole filled), con- 
sidering the ET dimers as unit cells the effective band- 
filling becomes ^. It has been suggested that pressure 
modifies the interdimer electron hoppings and changes 
the degree of frustration Experimentally, both 

an antiferromagnetic (AFM) state [18!| in systems with 
relatively large anisotropy, and a spin liquid state [l9j 
for nearly isotropic triangular lattices, have been found 
as well as SC. Theoretically, strong enough frustration 
within the Hubbard model destroys the AFM order, 
giving a paramagnetic metal (PM). Within theories of 
SC mediated by AFM fluctuations, the SC occurs in a 
narrow region sandwiched between the AFM and PM. 

0, i i Ea Ei E [S [3 m, m . 

We consider the Hamiltonian, 

H = -tY. + ^■^'■) - E + H.c.) 
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FIG. 1: (color online) Correlation functions as a function of 
U for t'=0.5. Solid line: bond order B' along the diagonal 
bonds. Dashed line: 5'CT(7r,7r). Both are discontinuous at the 
AFM-PM transition at U = 6.1. Lines are guides to the eye 
based on calculations performed on a grid of size Af/ = 0.1. 



In Eq.[T]c|^ creates an electron with spin a (t, |) on site 

*7 "n-i^tj = cl^Ci^a- U is the on-site Hubbard interaction. 
The lattice structure is the conventional square lattice 
(hopping integrals t) with additional hopping integrals t' 
across the x + y diagonals for a total of three bonds per 
site. The ratio t' /t is for the square lattice and 1 for the 
isotropic triangular lattice. In what follows, we express 
all quantities in units of t. 

Accurate solutions to Eq. [1] are difficult to reach us- 
ing existing analytical approaches. Many of the tech- 
niques that find SC are based on mean-field embedded 
cluster methods that use extremely small cluster sizes, 
including some as small as 4 or 8 sites [14, m. Eq. [1] 
and a variant with two diagonal hoppings have recently 
been investigated within a numerically more precise Path 
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Integral Renormalization Group (PIRG) technique [20|. 
The phase diagrams found by these authors contain PM, 
AFM, and nonmagnetic insulator (NMI) phases only and 
suggest implicitly the absence of SC along the AFM- 
PM boundary. The authors, however, did not calculate 
the SC pair-pair correlations. Recent exact diagonaliza- 
tion calculations also investigated the AFM-PM phase 
boundary, but did not calculate pair-pair correlations 
[2l| . Specifically because of the difficulty of calculating 
SC pair-pair correlations, and the prediction of SC occur- 
ring only over a narrow region of the phase diagram, it 
is essential that SC correlations are calculated explicitly 
with high precision over the complete phase space. 

We report precisely such calculations within Eq. [T] via 
exact diagonalization in the present Letter. We focus 
here solely on ^-filling, numerical convergence for which 
is reached much faster than in doped systems. Although 
exact diagonalizations cannot necessarily establish the 
occurrence of SC, they can test precisely the necessary 
conditions for SC. First, if the model does indeed have 
a SC ground state with pair-pair correlations exhibiting 
LRO in the thermodynamic limit, finite clusters should 
necessarily exhibit short range order (SRO) in the pair- 
pair correlations. Second, if SC driven by the Hubbard 
U interaction is present, then there must exist some pa- 
rameter region where U enhances the pair-pair correla- 
tions relative to U = 0. Note that for the one strongly- 
correlated model where it is known that SC exists, the 
negative-t7 Hubbard model, both the above conditions 
are met [g^]. We show here that neither of the above 
criteria are met within Eq. [l]for repulsive U at ^-filling. 
For all U and anisotropy, the SC pair correlations become 
monotonically weaker with increasing U . 

We solve exactly the ground state of Eq.[T]on a periodic 
4x4 lattice. We calculate the spin structure factor 5cr(q) 
and the bond order Bij between sites i and j, 

•^-(q) = ^E^''''^'""""'^«"Tr-^ir)(«Tr'-«ir')>(2) 
rr' 

B^3 = E(4^C,-,+^-C-) (3) 

In Eq. [2] is the number of lattice sites. We define the 
standard pair-creation operator A| 

^I = ;^E5H(4T4..i-4i4+.,T) (4) 

In Eq. [3]the phase factor g{y) determines the symmetry 
of the superconducting pairs. For s-wave pairing g(y) is 
+1 for the four nearest-neighbor sites i~\-x^ i + y, i — x and 
i — y; for d^2_y2 pairing g^u) alternates as +1,-1, +1,-1 for 
the same four sites; and for dxy pairing g{iy) alternates 
sign over the four sites i+x+y, i+x—y, i — x—y, i — x+y. 
We will calculate the pair-pair correlation function as a 
function of distance, P{r) = (A|Ai+f). 




FIG. 2: (color online) Correlation functions as a function of U 
for t'=l.l. (a) Diagonal bond order B' (solid line), Script, tt) 
(dashed line), and SaiT^jO) (dot-dashed line), (b) Derivative 
of B' and So-(vr, 0) with respect to U. Vertical lines show esti- 
mated phase boundaries between PM, NMI, and (7r,0) phases. 
Lines are guides to the eye and based on calculations per- 
formed on a grid of size AU = 0.1; derivatives calculated 
with a three-point centered difference formula. 



We first discuss the non-SC phases within the model. 
In the limit i' = 0, the ground state of Eq. [T]has AFM 
LRO for any nonzero U . This is usually determined from 
plots of 5*0- (q) versus q, which exhibit strong peaks at 
q = (tt, tt). For t' ^ 0, the ground state is AFM only 
when U exceeds a critical value Uc{t') and remains a 
PM for smaller U. We determine Uc from calculations of 
S'o.(vr,7r), the bond order corresponding to the diagonal 
t' bonds {B'), and the ground state expectation value 
of the double occupancy d — {ni^-^Ui^i). In Fig. [T] we 
plot ^^(Tr, tt) and B' for t' = 0.5. It is in this region of 
t' that the d^^-y^ SC phase is claimed to be broadest 
(l6| . Uc is identified from a discontinuous decrease in 
B' accompanied by a simultaneous increase in S'o-(7r,7r). 
The double occupancy d also decreases discontinuously 
at the same Uc (not shown). The large decrease of B' 
and d indicate loss of electron mobility associated with a 
metal-insulator (M-I) transition, while the sharp increase 
in S'(j(7r,7r) indicates the AFM nature of the insulator. 
The simultaneous discontinuity of all observables at Uc 
is consistent with a first-order (discontinuous) quantum 
phase transition. The discontinuity continues to exist for 
t' < 1. The size of the discontinuity in B' decreases as 
t' approaches 1; this is likely related to the appearance 
of the NMI phase as well as the change in the magnetic 
ordering at = 1 we discuss below. We also performed 
calculations (not shown here) on a 12 site cluster, chosen 
so that (tt, tt) AFM order is not frustrated at t' = 0, and 
found comparable values of Uc- 

The magnetic ordering is qualitatively different for 
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t' > 1 where Sa{<l) peaks at q = (tt, 0) instead of 
(tt, tt) [20|, with the amphtude of the peak increasing 
with U. This transition is expected from the known re- 
sults for the antiferromagnetic Heisenberg spin Hamil- 
tonian with nearest (Ji) and next-nearest (J2) neigh- 
bor exchange interactions (two diagonal bonds per site) 
[2^. Here J2/J1 = 0.5 is a quantum critical point, with 
q = (tt, tt) dominating for J2/J1 < 0.5, and q = (7r,0) 
for J2/J1 > 0.5 [IS]- This value of J2 corresponds to 
t' ~ within the U ^ 00 limit of the Hubbard 

model with two diagonal bonds, and this is where the 
transition to the (7r,0) ordering is found From per- 
turbation theory, the quantum critical J2/ Ji shifts to 1 
when the number of diagonal bonds per site is one in- 
stead of two, corresponding to t' = 1 within the large U 
Hubbard model. This is exactly where we have found the 
transition in our numerical work. 

In Fig. [2ta) we show B' and spin structure factor data 
for t' = 1.1 as a function of U. Unlike the discontinuous 
transition between the PM and AFM phases in Fig. [TJ 
the correlation functions change continuously now. Im- 
portantly, the steepest decrease in B' and the changes 
in S„{q), for both q — (7r,7r) and q — (tt, 0), occur at 
different U, suggesting the existence of a third phase be- 
tween the PM and the (tt, 0) magnetic order. In Fig.[2I^b) 
we show the derivatives of B' and ^^(Tr, 0) with respect 
to U. At U = 7.3 a minimum occurs in dB' /dU , which 
is accompanied by a minimum in the derivative of the 
double occupancy d at nearly the same value of U (not 
shown). At U = 13.8 similar maxima occur in derivatives 
of B' and S'CT(7r, 0). We therefore identify three distinct 
phases: a PM phase for U < 7.3, a nonmagnetic (NM) 
phase for 7.3 < U < 13.8, and the magnetically ordered 
(tt, 0) phase for U > 13.8. The transport behavior within 
the NM phase is difficult to determine from our finite- 
cluster calculations. The continuous decrease of B' and 
d in the intermediate U region is, however, consistent 
with an insulating character of this phase. We there- 
fore identify the intermediate U-phase as the previously 
discussed NMI phase 

Mean field calculations finding SC generally agree that 
the pairing is of d,j.2_y2 symmetry for Q < t' < 1. In 
Fig. [3^a) - (c) we show the dx2_y2 singlet pair-pair corre- 
lation function Pd ^ , (r) as a function of distance r and 
Hubbard interaction U, for t' — 0.2, 0.5 and 0.8, respec- 
tively. It is for the intermediate t' ~ 0.5 region that SC 
correlations have been claimed to be the strongest. The 
insets show the behavior of Pd 2 2 ^ ^ function of 

U for r — 2, and 2a/2. Pd 2_ 2 (^) the two larger 
values of r that we have chosen are least affected by the 
finite lattice size, since there is no overlap between pairs 
on the 4x4 lattice for these r. 

Several features in the pair-pair correlations are im- 
mediately apparent. First, the pair-pair correlations de- 
crease rapidly with distance, the values at the largest 




FIG. 3: (color online) (a)-(c) d^2_y2 pair-pair correlation 
functions as a function of r for t'—0.2, 0.5 and 0.8, respec- 
tively. Insets show [/-dependence of d^2_y2 pair-pair correla- 
tion functions for three different r. The PM-AFM transition 
is indicated by vertical lines, (d) U-dependence of the d^y 
pair-pair correlations at two different r for for t' = 1.1. Ver- 
tical lines mark PM-NMI and NMI-(7r, 0) transitions. Lines 
are guides to the eye and based on calculations performed on 
a grid of size AU = 0.1. 



r being more than an order of magnitude smaller than 
those at r = 0. Second, like the other correlation func- 
tions (Fig.[T]), Pd^2^ 2 ('') changes discontinuously at the 
PM-AFM transition. Finally, for r > 1, the pair-pair 
correlations are strongest at J7 = and decrease mono- 
tonically with increasing U . Pair-pair correlations at 
r — Q are trivially enhanced with increasing U . Nearest- 
neighbor correlations [r — 1) decrease with U for U < Uc 
for all t' but increase for U > Uc for large enough t'. 
This increase is associated with the transition to AFM 
and occurs only on the insulating side of the PM-AFM 
boundary. Furthermore, the small increase in the r = 1 
correlations is accompanied by a large simultaneous de- 
crease in the long-range component, consistent with the 
entry into an insulating state (clearly seen in the J7 = 4 
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FIG. 4: Phase diagram as a function of t' and U . Lines are 
guides to the eye. The PM-AFM boundary is located to 
within AC/ = 0.1 using the discontinuity of correlation func- 
tions. The PM-NMI and NMI-(7r, 0) phase boundaries are 
drawn from the inflection point in the correlation functions 
(see text.) 



data in Fig. EJa) and the f7 = 7 data in Fig. W^).) This 
behavior of the short-range correlations may be one rea- 
son embedded cluster methods using small cluster sizes 
find SC. 

Based on the q — (tt, 0) magnetically ordering for t' > 
1, it has been proposed that the SC correlations have 
pairing symmetry d^y or s -I- dxy [HI, [3, [111 ■ In Fig- E^d) 
we show the J7-dependence for dxy pair-pair correlations 
at t' — 1.1, for the two separations where the pairs are 
nonoverlapping. Exactly as for i' < 1, we find no region 
of U where the "long-range" d^y pair-pair correlations 

We have 
xy pairing 



are enhanced compared to their U — values 
performed similar calculations for s and s + d. 
symmetries, and find no enhancement for these either. 

We give the phase diagram from our calculations in 
Fig. 21 We have not plotted data for t' exactly 1 due to 
problems associated with degeneracy at this point. We 
have found that the NMI phase not only exists for t' > 1, 
but enters the phase diagram at approximately t' > 0.6; 
for all t' > 0.5 inflection points as in Fig. [2] appear in 
B' and d. Our phase diagram for the insulating states 
and the boundaries between them are in good agreement 
with other numerical results 1 IHIililia,^^]. We find, 
for example, that the minimum U necessary for the NMI 
phase is J7 ^ 7. This is very close to the boundary of f7 ^ 
6.7 at t' — 1 found recently using a variational cluster 
approximation [2^. The data in Fig. [H suggest that the 
PM-NMI transition is continuous [2u\, but confirmation 
of this would require finite-size scaling. 

Our most important result is that SC pair-pair correla- 
tions do not exhibit even SRO within Eq. 1, in any chan- 
nel and for any range of parameters. Theoretical models 
for SC in the CTS should either include electron-phonon 
interactions, or should be based on a different electronic 
model. We have recently suggested a correlated-electron 



theory of SC in the CTS that emphasizes their ;|-filled 
band nature, and the transition to the SC state from a 
density wave composed of local singlets [27| . Such a the- 
ory can explain naturally the recently observed transition 
from a valence-bond solid to the SC state @| 
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